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5 L'objet de I'invention conceme un procede de demodulation 

aveugle de signaux emis par plusieurs emetteurs et recus par un reseau 
compose d'au moins un capteur. 

II s'applique par exemple pour un reseau d'antennes dans un 

contexte electromagnetique. 
10 L'objet de I'invention conceme notamment la demodulation de 

signaux, c'est-a-dire I'extraction des symboles {a k } emis par un emetteur 
module lineairement. 

La figure 1 represente un systeme de traitement d'antennes. 

15 comportant plusieurs emetteurs El. el d'un systeme T de traitemeni 
d'antennes comportant plusieurs antennes Ri recevant des sources radio-: 
electriques avec des angles d'incidence differents. Les angles d'incidence| 
des sources ou emetteurs peuvent etre parametres soil en 1D avec l'azimu| 
9 m soit en 2D avec Tangle d'azimut e m et I'angle d'elevation A m . 

20 La figure 3 schematise un principe de modulation et demodulation 

des symboles { a k } emis par un emetteur. Le signal se propage au travers 
d'un canal a multi-trajets. L'emetteur emet le symbole a k a I'instant k.T, ou T 
est la periode symbole. La demodulation consiste a estimer et a detecter les 
symboles pour obtenir en sortie du demodulateur les symboles estimes a k . 

25 Sur cette figure, le train de symboles {a k } est filtre lineairement a remission 
par un filtre d'emission H appele egalement filtre de mise en forme h 0 (t). 

Dans la suite de la description, on definit sous I'expression 
« demodulation aveugle », des techniques qui n'utilisent pas d'information a 
priori sur le signal emis : filtre de mise en forme, sequence d'apprentissage, 

30 etc.. 
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Les dix dernieres annees ont vu le developpement des techniques 
de demodulation aveugle SIMO, abrege de entree unique sorties multiples, 
(en abrege anglo-saxon Single Input Multiple Output) dites a sous-espace 
utilisant les statistiques d'ordre 2, telles que decrites dans la reference [7]. 
5 Ces algorithmes pr^sentent toutefois Tinconvenient de ne pas etre robustes 
ni a une sous-estimation ni a une sur-estimation de I'ordre du canal de 
propagation : etalement tempore! dependant des multMrajets et du filtre de 
mise en forme. Pour contourner ce probleme, il a ete propose une technique 
de prediction lineaire decrite dans la reference [11] qui presente comme 

10 inconvenient d'etre moins performante lorsque la longeur du canal est 
connue. Pour ameliorer les techniques & sous-espace, la methode decrite 
dans [16] propose une technique parametrique qui necessite 
malheureusement la connaissance du filtre de mise en forme. 

Dans la reference [13], les auteurs proposent une technique a 

15 base de covariance matching, ayant notamment inconvenient d'etre tres 
difficile a mettre en ceuvre. C'est ainsi qu'il a ete developpe une technique 
sous-optimale decrite dans la reference [12] plus facile a mettre en oeuvre en 
minimisant un critere de vraisemblance et supposant le caractere gaussien 
des symboles.. Cette hypothese n'est pas verifiee pour les modulations 

20 lineaires couramment utilisees telles que les PSK (Phase Shift Keying) ou les 
QAM (abreviation anglo-saxonne de Quadrature Amplitude Modulation). 

II est aussi conhu dans les methodes CMA (en abrege anglo- 
saxon Constant Modulus Algorithm) d'utiliser une approche spatio-temporelle 
decrite par exemple dans la reference [6]. Cette famille de methodes 

25 presente toutefois I'inconvenient de n'etre adaptee qu'a une classe 
particuliere de modulations telles que les PSK qui sent a module constant, 
Cette methode est iterative et a done Tinconvenient de devoir etre 
correctement initialise. Pour finir, les methodes CMA presentent le 
desavantage de converger moins vite que la methode a sous-espace 

30 enoncee precedemment. D'autre part, la reference [20] decrit une methode a 
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sous-espace exploitant les statistiques d'ordres superieurs pour 
I'identification de canaux a Reponse Impulsionnelle Finie (FIR) et a phase 
non minimum. 

5 L'objet de la presente invention concerne un procede base 

notamment sur des techniques de separation de sources en aveugle 
connues de I'Homme du metier et decrites par exemple dans les references 
[4] [5] [15] [19] supposant que les symboles emis sont statistiquement 
independants. Pour cela le procede construit une observation spatio- 
10 temporelle dont les sources melangees sont des trains de symboles de 
I'emetteur. Chaque train de symboles est par exemple le meme train de 
symboles decale d'un nombre entier de periode symbole T. 

L'invention concerne un procede de demodulation aveugle d'une 
source ou emetteur de forme d'onde lineaire dans un systeme comportant 
15 une ou plusieurs sources et un reseau de capteurs et un canal de 
propagation caracterise en ce qu'il comporte au moins les etapes suivantes : 
. determiner le temps symbole T et on echantillonne a Te tel que T=ITe 
(I entier), 

. a partir des observations x(kTe), construire une observation spatio- 
20 temporelle z(t) dont les sources melangees sont des trains de 

symbole de I'emetteur , 
. appliquer une methode de type ICA sur le vecteur d'observation z(t) 
pour estimer les U trains de symboles { a m ., } associes aux vecteurs 

* A 

de canal \j=K( k j), 

25 • ordonner les U sorties (a m , h h„ ) dans le m§me ordre que les entrees 
(am-i , h 2 (i)) afin d'obtenir les vecteurs de canal de propagation 

• determiner la phase aimax associee aux sorties. 



Le procede selon l'invention offre notamment les avantages 

suivants : 
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• II ne fait aucune hypothese sur les constellations de symboles 
contrairement aux methodes decrites dans Tart anterieur, 

• II ne necessite pas la connaissance du filtre de mise en forme, 

• Le module des symboles n'est pas suppose constant, . 

5 • II est robuste a une surestimation de la longueur du canal, 

• II permet de traiter le cas des canaux de propagation avec des trajets 
correles, 

• II est direct et simple a mettre en ceuvre sans etape de recoupement 
des trajets correles. 

10 

D'autres caracteristiques et avantages de Pobjet de la presente 
invention apparaTtront mieux a la lecture de la description qui suit donnee a 
titre illustratif et nullement limitatif a la lecture des figures annexees qui 
represented : 
15 * La figure 1 un exemple d'architecture, 

• La figure 2 les angles d'incidence des sources, 

• • La figure 3 le processus de la modulation lineaire et demodulation 

d'un train de symboles, 

• La figure 4 le schema d'un emetteur a modulation lineaire, 

20 • La figure 5 un resume du principe general mis en ceuvre dans 
('invention, 

• La figure 6 la representation d'une constellation, 

• La figure 7 un premier exemple de mise en ceuvre du procede ou le 
signal est re<?u en bande de base, 

25 * La figure 8 un deuxieme exemple ou le signal est re?u en bande de 
base et les multi-trajets sont d^correles, 

• La figure 9 un troisieme exemple oCj le signal est re^u en bande de 
base et les multi-trajets sont decorretes par groupe. 
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Afin de mieux faire comprendre le procede selon I'invention, la 
description qui suit concerne un procede de demodulation aveugle aux 
ordres superieurs d'un Smetteur de forme d'onde lineaire dans un reseau 
ayant une structure telle que celle decrite a la figure 1 , par exemple. 

5 Avant d'expliciter les etapes mises en ceuvre par le procede\ on 

decrit le modele du signal utilise 
Modele du signal emis par une source ou emetteur 
Modulation lineaire 

Les figures 3 et 4 montrent le processus de la modulation lineaire 

10 d'un train de symboles {a k } a la cadence T par un filtre de mise en forme 
h 0 (t). 

Le peigne de symboles c(t) est tout d'abord filtre par le filtre de 
mise en forme h 0 (t) et ensuite transpose a la frequence porteuse f 0 . Le filtre 
NRZ, qui est une fenetre temporelle de longueur T , tres souvent defini par 

15 ho(t)=nT(t-T/2), est un exemple particulier non limitatif de filtre Remission. 
Dans les radiocommunications, il est egalement possible d'utiliser le filtre de 
Nyquist dont la transformee de Fourier h 0 (f)=n B (f-B/2) se rapproche d'une 
fenetre de bande B, lorsque le roll-off est nul alors h 0 (f)=n B (f-B/2) (le roll-off 
d^finit la pente du filtre en dehors de la bande B). 

20 Le signal module s 0 (t), emis par I'emetteur, s'ecrit a I'instant t k =kT 8 

(T e : periode d'echantillonnage) en fonction du peigne de symboles c(t) : 

(1). 

s 0 (kT e )=X ho(iT e ) c((k-i)T e ) 

Prenons un temps symbole T egal a un nombre entier de fois la 
periode d'echantillonnage, T=IT e et posons k=ml+j avec 0 < j < I. Puisque 
c(t)=Z r a r 5(t-rlT e ), autrement dit, comme c(t)=a u pour t=ulT e et c(t) = 0 pour 
25 t*ulT e , les seules valeurs de i pour lesquelles c((k-i)T e ) est non nul veYifient k- 
i=ul, c'est-a-dire telles que i=ml+j-ul=nl+j ou n=m-u. Finalement, I'expression 
(1) devient : 
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& (2) 
s 0 (mlT e + jT e )= 2] h 0 (nlT e +jT e ) a m . n pour 0 < j < I 

Le parametre L 0 est la demi-longueur du filtre d'emission qui s'etale sur une 
duree de (2L 0 +1)IT e . Dans le cas particulier d'un filtre d'emission NRZ, on 
obtient L 0 =0. Quant au signal emis s(t) f il verifie s(t)=s 0 (t) exp(j27tf 0 t) car il est 
egal au signal s 0 (t) transpose a la frequence f 0 . Dans ces conditions 
5 ('expression de s(mlT e +jT e ) est d'apres (2) : 

s(mlT e +jT e )= £ ho(nlT e +jT e ) expG27if 0 (nl+j)Te) exp(j27rf 0 (m-n)iT e ) ^ 

«=-/„ 

u 

= £ h F0 (n I T e + jT e ) be™ tel que 0 < j < I 
ou h F o(iT e )=h(iTe) exp(j27tf 0 iT e ) et bi=aj exp(j27if 0 ilT e ) 



Reception des signaux sur les capteurs 

Le signal s(t) (FIG.3) emis passe k travers un canal de propagation 
avant d'etre re<?u sur un reseau compose de N antennes. Le canal de 
10 propagation peut se modeliser par P multi-trajets d'incidence e p , de retard t p 
et d'amplitude p p (1 < p < P). En sortie des antennes on a le vecteur x(t) qui 
correspond a la somme d'un melange lineaire de P multi-trajets et d'un bruit 
suppose blanc et gaussien. Ce vecteur de dimension M<1 a I'expression 
suivante : 

p (4) 
x(t) p P a(9 p ) s(t-T p )+ b(t)=A s(t) + b(t) ' ' 

15 ou p p est Pamplitude du p ldm9 trajet, b(t) est le vecteur bruit suppose 
gaussien, a(0) est la reponse du reseau de capteurs a une source 
d'incidence 9, A=[ a(9i)... a(0p)] et s^Hs^-t^^s^-tp)] 7 . En notant que x p = 
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r p T+Axp (ou (0<AT p <T=IT e ) et r p est un entier) et en utilisant I'expression (3) 
dans I'equation (4), on obtient pour le vecteur sur les antennes : 

x(mlTe+jT e ) =2 2 P P a(6 p ) h F0 (n IT e +jT e -AT p ) b m . n . rp + b(mlT e +jT e ) 

En effectuant le changement de variable suivant u p = n + r p , le vecteur regu 
par les antennes s'exprime : 

p (6) 
x(mlT e +jTe) =]T 2 Pp a ( Q p) h Fo((u P - r p )ITe+jT e -ATp) b m _ Up + b(mlT e +jTe) 

5 En notant r min = min{r p } et w = max{r p }, I'equation (6) peut etre reecrite de la 
maniere suivante : 

x(mlT e +jT e ) = £ ^ P P a(6 p ) h F0 ((u- r p )IT e +jT e -AT: p ) lnd [rp . L o,rp + LO](u) b m _ u ( * 
+ b(mlT e +jT e ) 

Ou lnd [r ,q](u) est la fonction indicatrice usuelle (lnd [r , q ](u)=1 pour r<u<p et ; . .. 
Ind [f>q] (u)=0 autrement) definie sur I'ensemble des entiers relatifs a valeur ^ 
dans I'ensemble binaire {0, 1}, caracterisee par lnd lr , ql (u) = 1 si u appartient a " 
10 I'intervalle [r,q] et lnd (r , ql (u) = 0 sinon. De ce fait, en notant v(t) le vecteur de 
canal : 

p (8) . 

v(ulT e +jT e ) =2 PP a ( 0 p) h F0((U" r p )IT e +jT e -AT p ) lnd trp -L0. rp+L0](U) 

p=\ 

ou t= ulT e +jTe et I'expression (5) devient : 

w*. (9) . 

x(mlTe+jT e )= £ v(ulT e +jT e ) b m . u + b(mlT e +jT e ) 

Interference entre symboles 

Le vecteur observation x(t) issu du reseau d'antennes a I'instant 
15 t = mlTe+jT e fait, d'apres I'equation (9), intervenir le symbole b m mais 
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egalement les symboles b m - u ou u est un entier relatif appartenant a 
Tintervalle [r min -Lo, r ma x+Lo], phenomene qui est plus connu sous le nom 
d'lnterference Entre Symboles (IES). Notons L c le nombre de symboles 
participant a TIES et bornons Tintervalle de valeurs prises par ce dernier. 
5 D'apres ('equation (9), si Intersection des intervailes [r p -Lo, r p +L 0 ] est non 
vide, alors on a L c = | r max - r m j n l + 2L 0 + 1 . De ce fait, lorsque r max = r m j n , 
c'est-a-dire lorsque, tous les multi-trajets sont correles, la borne minorante de 
U est atteinte et vaut L c = 2L 0 + 1. Ce cas se traduit egalement 

mathematiquement par) max{x p }-min {x p } | < T. Par contre, si ^intersection 

p p 

10 des dits intervailes est vide, et que le cas echeant, tous les intervailes [r p -L 0 , 
r P +U] sont disjoints, alors on a Lc==Px(2L 0 +1), ce qui constitue une borne 
majorante a ('ensemble de valeurs susceptibles d'etre prises par L c . Ce 
dernier cas de figure correspond concretement au cas de multi-trajets tous 
decorreles deux a deux, ce qui mathematiquement peut egalement s'ecrire 

15 " Vi 9ft j, | n — rj |>2Lo , condition obtenue des que | x\ - Tj | > (2L 0 + 1)T. Pour 
resumer, la quantite L c verifie de maniere generate Tencadrement suivant : 



Uexpression traduisant le vecteur re?u par les capteurs peut alors se reecrire 
de la maniere suivante, ou cette fois n'apparaissent que les L c symboles b m .y 
d'interet : 



2L 0 +1 < U< Px(2L 0 +1) 



(10) 



X(mlT e +jT e ) =2J h(n(l)IT e +]T e ) b m . n (i) + b(mlT e +jT e ) 



(11) 



20 Ou V 1 <I<U, et rmin-L 0 < n(l) < r min +L 0 et ou: 



h(t)«£ p p a(e p )h F0 (t-T P ) 



(12) 



Techniques ICA 

Le procede fait appel a des techniques ICA basees sur le modele 
suivant donne a titre iliustratif et nullement limitatif : 
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* ^ (13) 

Uk=2_, 9 j s ik + n k = G s k + n k 
/=i 

ou u k est un vecteur de dimension Mx1 recu a I'instant k, s ik est la i' eme 
composante du signal s k a I'instant k, n k est le vecteur de bruit et G =[ gi ... 
gj. Les methodes ICA ont pour objectifs d'extraire les l=L composantes s ik et 
d'identifier leurs signatures g-, (La reponse vectorielle de la source i au travers 

5 de I'observation u k ) a parti r des observations u k . Le nombre l=L de 
composantes doit etre inferieur ou egal a la dimension M du vecteur 
d'observation. Les methodes des references [4] [5] et [15] utilisent les 
statistiques d'ordre 2 et 4 des observations u k . La premiere etape utilise les 
statistiques d'ordre 2 des observations u k (ces observations peuvent etre 

10 fonctions des signaux recus sur les capteurs) pour obtenir une nouvelle 
observation Zk telle que : 

l „ (14) 

z k = WiUk=2] gi Sik +n k =G s k + n k 
1=1 

ou les signatures & (1<i<L) sont orthogonales, G =[gi ... gt] et s k =[ s 1k ... 
s Lk ] T . La deuxieme etape consiste a identifier la base orthogonale des G a 
partir des statistiques d'ordre 4 des observations blanchies z k . Dans ces 
15 conditions on peut extraire les signaux s k en effectuant : 

s k = G* z k = G # Wi u k (15) 
Ou s k est I'estimee des signaux s k et ou # est I'operateur de pseudo-inversion 
definiparG* = (G H G)' , G H . 

La methode ICAR [19] utilise quant a elle uniquement les statistiques d'ordre 
4 pour identifier la matrice G = [gi ... g«] des signatures. 
20 En resume, I'idee mise en ceuvre dans le procede selon I'invention 

est de construire une observation spatio-temporelle dont les sources 
melangees sont des trains de symboles de I'emetteur. Chaque train de 
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symboles est par exemple le meme train de symboles decale d'un nombre 
entier de periode symbole T. 

Le procede decrit ci-apres comporte plusieurs variantes de 
realisation dont certaines sont expliquees a titre illustratif et nullement 
5 limitatif. 

Premiere variante de realisation du procede 

La figure 7 represente un premier exemple de variante de 
realisation du procede ou le signal est regu en bande de base. 

Le procede comporte une etape 1.1 de determination du temps 
10 symbole Te en appliquant par exemple un aigorithme de detection cyclique, 
tel que celui decrit par exemple dans [1] [10]. 

L'etape suivante I.2 consiste k interpoler les observations x(t) a / 
echantillons par symbole, tel que T=fTe. 

Dans ces conditions ou f 0 =0 et b k =a k , Pexpression (11) du vecteur devient : 

l< (16) 
x(ml T e + jT e ) =2 h(n(l) I T e + JT e ) a m - n (i) + b(ml T e + jT e ) pour 0 < j < I 

Comme Tequation (16) est verifiee pour 0 < j < /, le procede construit 
Tobservation spatio-temporelle suivante (etape I.3) a partir des observations 
x(kTe): 



15 



z(ml T e )= 



x(mIT e 



=2 h z (n(l)) a m .n(i) + b z (ml T e ) ou h 2 (n)= 



(17) 



avec h nj = h(nl T e + jT e ) et b 2 (ml T e )=[ b(ml T e ) T ... b(ml T e +( l-1)T e ) T ] T - 
Sachant que x(t) est de dimension Nx1 , le vecteur z{t) est alors de dimension 
20 Nlx1. 

h(k) est un vecteur dont la n i6me composante est le k ,6ma coefficient du filtre 
filtrant lineairement le train de symbole {a m } sur le n ,Sme capteur. Le filtre de 
coefficient vectoriel h(k) depend a la fois du filtre de mise en forme et du 
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canal de propagation. 

Afin d'extraire les U trains de symboles { a™ } d'interet (nombre 
de symboles qui participent a I'lES), le procede echantillonne le signal recu a 
l=(2L 0 +1), en supposant que P < N. 
5 Sachant que le filtre NRZ verifie 2L 0 +1=1 et le filtre de Nyquist 

2Lo+1=3 pour un roll-off de 0.25, les trains de symboles peuvent etre extraits 
pour ces deux filtres de mise en forme respectivement lorsque P < Nl et 3P 
<NI . 

Le vecteur d'observations z(t) etant determine, le procede 
10 applique une methode de type ICA pour estimer les U trains de symboles 
{ a m -i } associes aux vecteurs de canal h sj =h.(kj). 

La p me sortie des methodes ICA donne le train de symboles {am.j } 
associe au vecteur de canal h,,, . Les trains de symboles {a m . j } estimes 
arrivent dans un ordre different de celui des trains { a„vi } en verifiant : 

a m , j = p expGcci) a m -i et h 2j = h z (i) C 8 ) 

15 Les trains de symboles { a m , j } sont estimes avec la meme amplitude car les 
trains de symboles {a m --,} sont tous de m§me puissance en verifiant : 

E[|am-n(1)| 2 ]=---=E[| a fn-n(Lc)| 2 ]- 

L'etape suivante I.4 du procede a pour objectif d'ordonner les U 
sorties (am, j, h..,,) dans le m§me ordre que les entrees (a m .i , h 2 (i)) afin 
20 d'obtenir les vecteurs de canal h,.,=h,(A y ) . Pour cela, le procede intercorrele 
deux a deux les sorties a mii et ^ en calculant le critere c it j(k) suivant : 
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Lorsque la fonction |cij(k)| est maximum en k=k max les i i6me et j i§me sorties 
verifient ; a m ,i= a m .kmaxj. L'algorithme de classement des sorties a m , n (i) 
a m ,n(Lc> est par exemple compose des etapes suivantes : 
Etape n°A.1 : Determination de la sortie a m ,i max associee au vecteur de canal 

5 h :Jmax de plus fort module. 

Etape n°A.2 : Pour toutes les sorties & m -Ki ou ]Vi max determination des indices 
k=kj maximisant le critere | Ci ma x t j(k)j. On en deduit pour chaque j que a m , {max = 
a m - k jj. Sachant que c imax> j(k|)= exp(jotj max - jcq) on remet la j i6me sortie a la meme 
phase que la i m ax idme sortie en effectuant a m - kj = c imaX| j(kj) a m , On remet aussi 

10 en phase les vecteurs de canal en effectuant : h r (fc,)= h ZJ Cj maXtj (kj)\ 

Etape n°A.3 : Cette etape remet en ordre les sorties a m . kj et les vecteurs de 
canal h :J =h s (kj) dans Tordre croissant des kj sachant que a m = S m jmax et que 

h ? (0)=h Ztimax . 

En sortie de ces trois etapes, on obtient les trains de symboles 
15 { a m - k } associts aux vecteurs de canal h,(* y )' Sachant que les symboles 
estimes verifient a m . k = exp(j a i max ) a m . k , la demfere 3tape du proctde va 
consister a estimer cette phase ccj max . Pour cela on identifie tout d'abord la 
constellation des symboles a k parmi une base de donnees composee de 
Pensemble des constellations possibles. Cette base est constitute des 
20 constellations connues tel que nPSK, n-QAM. A chaque fois que Ton 
detectera ou que Ton aura connaissance d'une nouvelle constellation on 
enrichira la base. 

La figure 6 represente un exemple de constellation de 8-QAM 
lorsque aj max =0 et ajm ax *0. Dans cet exemple de mise en ceuvre, le procede 
25 comporte alors les etapes suivantes : 

L'6tape 1. 5 suivante consiste a determiner la phase de sortie 
assoctee au vecteur de canal de plus fort module. Pour identifier la 
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constellation et determiner la phase le procede effectue par exemple les 

etapes suivantes : 

Etape 1 .5 = Etapes B.1 , B.2 et B.3 

Etape n°B.1 : Estimation des positions des etats de la constellation (points 
5 rouges sur la figure) par la recherche des maximums de I'histogramme 2D 
des points M k =( reel(a k ), imag(a k ) ). Pour une constellation a M etats on 
obtient M couples (Qm, v J pour 1<m<M. 

Etape n°B.2 : Determination du type de la constellation en comparant la 
position des etats (u m ,vj de la constellation des {a k } a une base de donnees 
10 composee de ('ensemble des constellations possibles. La constellation la 
plus proche est composee des etats (u m ,v m ) pour 1<m<M. 
Etape n°B.3 : Determination de la phase cw en minimisant au sens des 
moindres carres le systeme d'equations suivant : 

u m = cos(Oi m ax) u m - sin(cw) v m et v,, = sin(cw) u m + cos(cw) v m 

15 pour 1<m<M 

Le procede peut comporter une etape d'estimation des 

parametres du canal de propagation en angle 9 P et retard x p de I'equation £8) 
par I'algorithme propose dans [8]. L'etape consiste a extraire dans un 
premier temps les vecteurs h(nl T e + jT e ) pour 0<j<l des vecteurs de canal 
20 des vecteurs h>,) defini dans I'equation (17). Puis de construire la 
matrice H=[ h(n(1)l T e )... h(n(U) I T e )] de (11) avec les h(nl T e + jT e ) pour 
appliquer la methode [8] d'estimation parametrique des multi-trajets : (6 P , x p ) 
1<p<P. 

Deuxieme variante de realisation du procede 

25 Les figures 8 et 9 schematised une autre variante de realisation 

pouvant comporter deux variantes correspondant respectivement au cas des 
multi-trajets decorreles et au cas des multi-trajets correles par groupe. 
Cas des multi-trajets decorreles. 
Le signal est recu en bande de base avec { b k }={ a k } 
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Les multi-trajets dont les retards verifient | tj-ti | > (2L 0 +1) T, ont 
ravantage d'etre decorreles entre eux en verifiant : E[s(t-Tj) s{t-rj)*]=0. En 
observant l'equation (4), on constate alors qu'il suffit d'appliquer une 
methode de type ICA lorsque P<N sur I'observation x(t) pour obtenir les 
signaux s(t-x p ) de chacun des multi-trajets. Apres ('estimation des signaux 
des differents multi-trajets, le procede determine leurs puissances pour 
garder le signal s(t-T pma x) du multi-trajet de plus forte amplitude p pma x. Pour 
determiner ce trajet principal on utilise le fait que les sorties des methodes 
ICA verifient en asymptotique : 

(20) . 



x(t)=£ p p a(9 p ) s(t-T p ) = £ ft«i,(0 avec 5,(0= j(f *' } et 
&\=4Yp P? a(0 p ) 

10 ou Yp=p p 2 E[|s(t-x p )J 2 ]. Comme les vecteurs a(8 p ) sont norrnes en verifiant 
a(0 p ) H a(9 p )=N f le trajet d'amplitude maximum sera associe a la i m ax feme sortie 
OU OCimax^ Simax ^imax est maximum. Comme d'apres l'equation (3) la sortie 
W(')=s(t-T pma x) verifie : 

^(m IT e +jT e )= £ h F0 (n I T e + jT e -T pmax ) a m . n telque 0<j<l (21) 



on peut constituer le vecteur observation suivant : 



z(ml T e )= 



ls imax (mIT e + (J-l)T e ) 



= £ h z (n) a m . n ou h z (n)= 



h 



'n.l 



(22) 



15 ou h n ,j=h F o(n I T e + jT e -T pmax ). D'apres le modele de l'equation (22), il suffit 
d'appliquer une methode de type ICA sur I'observation z(ml T e ) pour estimer 
les 2U+1 trains de symboles { a,™ } avec - Lo<n<Lo- Pour extraire les 
incidences 0 P du canal de propagation, il suffit d'apres (20) de chercher pour 
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chaque signature ai (1<i<P) le maximum du critere c(9)=| a(9) H aj | 2 . Pour 
extraire les retards %■> -ti du canal de propagation, il suffit d'apres (20) de 
chercher pour chaque signal s s (t) (1<i<P) le maximum du critere 

c(xHW-*)Wl l - 

5 En resume cette variante comporte par exemple les etapes 

suivantes : 

Etape n°ll.a.1 : Determination du temps symbole T en appliquant un 
algorithme de detection cyclique comme dans [1][10]. 

Etape n°ll.a.2 : Echantillonnage des observations x(t) a / echantillons par 
1 o symbole tel que T=/ T e . 

Etape n°ll.a.3 : Application d'une methode ICA sur les observations x(t) pour 

obtenir s,(0 et aj pour 1<i<P. 

Etape n°ll.a.4 : Determination de la sortie i= imax ou ctj= a" a t maximum. 
Etape n°ll.a.5 : Constitution du vecteur observation z(t) de (22) a partir du 
15 signal s lmax (t). 

Etape n°ll.a.6 : Application d'une methode ICA pour estimer les trains de 
symboles { a m . n } oti -U<n< U- On choisit parmi les trains de symboles celui 
qui est associe au vecteur h 2 (n) de plus fort module : {a m }. 
Etape n°ll.a.7 : Determination de la phase cw de la sortie associee au 
20 vecteur h z (n) de. plus fort module en appliquant les etapes B.1 , B.2 et B.3. 
Etape n°ll.a.8 : Remise en phase du train de symboles {a m } en effectuant 
a m = a m exp(-ja im ax) . Le train de symboles {dj constitue la sortie du 
demodulateur de ce sous-procede. 

Etape n°ll.a.9 : Estimation des parametres du canal de propagation en angle 
25 9 P et retard t p en maximisant pour 1<i<P les criteres | a(G) H a\ | 2 et 
[^(/-r) | 2 pour respectivement les angles et les retards. 
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Cas general avec des multi-trajets quelconques ou correles par groupe 

Dans cette variante dont le schema est donne en figure 9, le 
procede considere qu'une partie des multi-trajets sont correles. En 
considerant que I'emetteur est re9u suivant Q groupes de multi-trajets 
correles, le vecteur signal recu par les capteurs de I'equation (4) devient : 

(23) 

Aq £2 q S(t, Tq) +D(t) 



X W =f] PP.q a ( G P.q) S(t-Tp.q)+b(t)= £ Aq Q q S(t, Tq) + 



OCl A q =[ a(9 1>q )... a(e P q,q)], Qq=diag([p 1iq ...p P q,q]) et S(t, Tq)=[ S(t-T,,q)... s(t- 

Tpq,q)] T avec s q =[ti,q ...T P q.q] T . En appliquant une m^thode ICA on estime en 
sortie du separateur les signaux et les signatures suivantes d'apres : 

(24) 

• 1 - X- ✓ -ivy 

A=[a!... a P Q, Q H A, U n ... A Q U Q ] U et s(/)=n 







" *,(0 " 


y Q s(m q )_ 







ou n est une matrice de permutation, U q V q =£2 q et V q E[s(t, x q ) s(t, x q ) H ] V q H 
=Ip q . Ainsi les trajets decorreles tel que E[s(t-x p , q ) s(t-x p * jq ')*]=0 sont regus sur 
des voies 5 f .(r) et 5 y (0 differentes. Les trajets correles ou E[s(t-x p , q ) s(t- 

x P ', q )*]*0 sont melanges sur une m§me voie et sont presents sur P Q k la 

fois. Dans !a 1 i§re etape de ce sous-procede nous utilisons ce resultat pour 
identifier les Q groupes de multi-trajets correles. En prenant les sorties i et j 
15 du separateur, les deux hypotheses suivantes peuvent §tre testees : 



Ho 



*(0 = *,(') 



f4(0 = a,5(r-r p ) + ^(0 

(0 



(25) 



ou E[bj(t) bj(t-T)*]=0 quelque soit la valeur de x. Ainsi dans I'hypothese H 0 il 
n'existe pas de multi-trajets communs aux deux sorties i et j et dans 
I'hypothese Hi il y en a au moins un. Le test va consister a determiner si les 
sorties £,(/) et f,(/-r) sont correlees pour au moins une des valeurs de x 

20 v^rifiant | x | < x max . Pour cela on peut appliquer le test de Gardner [3] qui 



17 



compare a un seuil le rapport de vraisemblance suivant : 

| ,2 K (26) 

V ii (x)=-2Kln(1- ) avec /v(*)=il/. (<)*,('-')' 

Ou Vij(x)<r| => hypothese H 0 
EtVij(x)>Ti => hypothese 

Le seuil T| est determine dans [3] par rapport a une loi du chi-2 a 2 degres de 
liberie. On cherchera tout d'abord les sorties associees avec la 1 ,ere sortie en 
lancant le test pour 2<j<P Q xQ et i=1 . Puis de la liste des sorties on enlevera 
5 toutes celles associees a la 1 Ifere qui constituera le 1 ier groupe avec q=1. On 
recommencera la meme serie de tests avec les autres sorties non correlees 
avec la 1 i6re sortie pour constituer le 2 i6me groupe. On effectuera cette 
operation jusqu'au dernier groupe ou au final il ne restera plus aucune voie 
de sortie. On obtiendra finalement en sortie du tri : 

Aq=A q U q et s 9 (0=V q s(t, x q ) pour (1<q<Q) (27) 

10 Les incidences 0 p ,q sont determinees a partir des A q pour (1<q<Q) 

en appliquant I'algorithme MUSIC [1] sur la matrice A q A q H . De ces 
goniometries on deduit les matrices A q . Sachant que. x q (t)= A q s,(/)= A q Q q 
s(t, x q ), on en deduit s(t, x q ) a une matrice diagonale pres en effectuant 
s(t,iq)= Aq%(t). Comme les elements des s(t,x q ) sont composes des signaux 

15 s(t-x p. q ), on determine les retards Xp.q-xi.-, en maximisant les criteres 
c(x)=|^. P (^-^) *u(0* I 2 ou est la P i6me composante de s(t,x q ). 
Sachant que E[s,(/) s 9 (0 H ]=lp q . que A q H A q = N l Pq et que A q s,(r) = A q a q s(t, 
x q ), on en deduit que le groupe de multi-trajets associes aux amplitudes £2q 
les plus importantes maximise le critere suivant : cri(q)=trace(A q A q ). On en 

20 deduit alors la meilleure sortie associee a A qm ax et s, mas (/). Comme d'apres 
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('equation (3) le vecteursCt, Iqmax ) verifie : 

*(/»/r tf+y x-r vmaXil ) 



s(mlT e+ jT e ,T £|max ) = 



s(m/T + ;T, - r ^ 



u (28) 
- Z h Fo(n'T e +jT eiIqmax )a m . n 



pour 0<j<letavec h F0 (n I T e+j T e , T qmax )= 



on peut constituer le vecteur observation suivant d'apres (27): 



z(mJ T e )= 



s qmax (mIT e + TJ 

qmax 



= Z h 2 (n)a m . n ou h 2 (n)= 



Vol 
h ' 



(29) 



ou h nj ._ Vqmax hFo(n , Te+ jTeIqmax) D , apr , s |e mo ^ ie ^ 

suff,t d'appHquer une metnode de type ICA sur ."observation z(m. T e ) pour 
5 est.mer les 2L 0+ 1 trains de symbole { a m . n } tel que - L 0 <n<L 0 . 

En resume cette variante comporte fes etapes suivantes • 

TV°T : D ' termlnati ° n du te ^ ™e T en appiiquant un 
algonthme de detection cyciique comme dans [1]fl 0 ] 

10 s^ro (e 7; b ' 2 : T EChantf " 0nn ^ des Nations x(t) a , echanti.ions par 
io symbole tel que T=l T e . 

Etape n-n.b.3 ; Applioation tfune ^ ^ ^ 

obtenir s(i) et A de I'equatiori (24). 

oTS' b t : T " des sorties suivant 0 9roupes de ™ m <*>° 

pourob.en.rA, e, s , (0 pour (1 £qsQ) ; Pour cela te8 , de de ^ 

.5 les couples de sodies i e. j par .e ,es, a deux hypotheses de PequaHon (26, 
On cherchera ,ou1 d'abord ,es sodies assoolees aveo ,a 1- sodie en JZ 
le .es, p 0ur 2< j2 P QX Q Bt w . Puis de , a (fete des sort . es 
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celles associees a la 1 fere qui constituera le 1 ier groupe avec q=1. On 
recommencera la m§me serie de tests avec les autres sorties non correlees 
avec la 1 i§re sortie pour constituer le 2 l6me groupe. On effectuera cette 
operation jusqu'au dernier groupe ou en final il ne restera plus aucune voie 
5 de sortie. 

Etape n°ll.b.5 : Determination du meilleur groupe de multi-trajets ou 
trace(A q H A q ) est maximum en q= qmax. 

Etape n°H.b.6 : Constitution du vecteur observation z(t) de (29) a partir du 
signal s qmzx (t). 

10 Etape n°H.b.7 : Application d'une methode ICA pour estimer les trains de 
symboles {a m . n } ou -L 0 <n< L 0 . On choisit parmi les trains de symboles celui 
qui est associe au vecteur h z (i) de plus fort module : {a m -i}. 
Etape n°ll.b.8 : Determination de la phase ow de la sortie associee au 
vecteur h 2 (i) de plus fort module en appliquant les etapes B.1 , B.2 et B.3. 

15. Etape n°H.b.9 : Remise en phase du train de symboles {a m } en effectuant 
l m = a m exp(-ja ima x) . Le train de symboles {<§„,} constitue la sortie du 
demodulateur de ce sous-procede. 

Etape n°H.b.10 : Estimation des parametres du canal de propagation en 
angle 9 q>p et retard T q>p . Les incidences e p , q sont determinees a partir des Aq 
20 pour (1<q<Q) en appliquant Talgorithme MUSIC[1] sur la matrice A q A q H . De 
ces goniometries on en deduit les matrices Aq pour en deduire une estime de 
s(t, x q ) en effectuant s(t,T q )= A q # x q (t). Comme les elements des s(t,T q ) sont 
composes des signaux s(t-T p .q), on determine les retards T P ,q-Ti,i en 
maximisant les criteres c(x)=\s qiP (t-r) 5,,,(0* | 2 ou s q ^t) est la p i6me 

25 composante de s (t,Tq). 

Autre variante de mise en ceuvre du procede 

Estimation de la frequence porteuse et deduction des {a m }. 

Cette technique va consister a estimer la frequence porteuse f 0 de 
I'emetteur ou le complexe z 0 = exp027if o T e ) pour ensuite dSduire les symboles 
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{a m } des symboles {b m } en effectuant d'apres (3): 

a m =b m exp(-j27if 0 mlT e )= b m z 0 " ml (30) 

Cette etape est appliquee apres I'etape n°l.4 de remise en ordre des 
symboles et des vecteurs de canal. D'apres les equations (3)(17) et (7)(8), 
on dispose des vecteurs de canal suivant : 

(31) 



h.(«) 



z* h(«/r c ) 
0 n '* l h(nIT c +T e ) 



L V' +(/ ~°h("/7; + (7-1)7;) 



pour ne£2 



ou Q. ={lnd( rp .uo, r P +L0](n)=1 pour un p tel que 1<p<P} 

5 Sachant que Q. ={ ni <...<nL C }. on constitue a partir des vecteurs h,(«)un 
grand vecteur b tel que : 



w = 



KM 



(32) 



.KM. 

La recherche de f 0 va consister a maximiser le critere suivant 
Porteuse(f 0 )=|w H c(expG27rf 0 T e ))| 2 



(33) 



Avec c(zq)= 



c(»i»^o) 
c(" 2 ,Zo) 

c(%»^o) 



et ou c(n,2 0 ) = 



!»/+(/-!) 
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Les etapes du procede adapte au cas d'un emetteur a frequence non nulle 
sont les suivantes : 

Etape n°lll.a.1 : Etapes 1.1 jusqu'& I.4 decrites ci-dessus pour obtenir les 
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trains de symboles {b m . kj } associes aux vecteurs de canal »...(*,) • 

Etape n°lll.a.2 : Construction du vecteur w de I'equation (32) a partir des 

Etape n°IILa.3 : Maximisation du critere Porteuse(f 0 ) de i'equation (33) pour 
5 obtenir fo- 

Etape n°JH.a.4 : Application de I'equation (30) pour deduire les symboles {a m } 
des symboles {b m }. 

Etape n°lll.a.5 : Etapes I.5 jusqu'a I.7 precedemment decrites. 

Dans le cas d'un emetteur a frequence non nulle et pour des multi- 
10 trajets decorreles les etapes sont les suivantes : 

Etape n°III.b.1 : Etapes ll.a.1 jusqu'a ll.a.4 decrites ci-dessus pour obtenir le 
vecteur z(t) de I'equation (22). 

Etape n°lli.b.2 : Application des methodes ICA [4] [5] [15] [19] pour estimer 
les Lc trains de symboles {b mJ } associes aux vecteurs de canal h„ . 
15 Etape n°IH.b.3 : Remise en ordre des trains de symboles {b mJ } et des 
vecteurs de canal fi^ en appliquant les etapes A.1. A.2 et A.3 pour obtenir 
les trains de symboles {&„_*, ) associes aux vecteurs de canal h rJ =h..(* ; ) • 
Etape n°lll.b.4 : Construction du vecteur w de ('equation (32) a partir des 

20 Etape n°lll.b.5 : Maximisation du critere Porteuse(f 0 ) de I'equation (33) pour 
obtenir fo. 

Etape n-lll.b.6 : Application de ('equation (30) pour deduire les symboles 
{a m } des symboles {b m }. 

Etape n°lll.b.7 : Choix du train de symboles associe au vecteur h z (i) de plus 

25 fort module : {a m -i}. 

Etape n°IU.b.8 : Etapes ll.a.7 jusqu'a ll.a.9 decrites precedemment. 

Dans le cas d'un emetteur a frequence non nulle et pour des multi-trajets 

correles les etapes sont par exemple les suivantes : 
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Etape n°III.c.1 : Etapes Il.b.1 jusqu'a ll.b.6 n°2.2 pour obtenir le vecteur z(t) 
de ('equation (29). 

Etape n°ll!.c.2 : Application des methodes ICA [4] [5] [15] [19] pour estimer 
les L c trains de symboles {b m j } associes aux vecteurs de canal h : j . 

5 Etape n°llf.c.3 : Remise en ordre des trains de symboles {b m J } et des 
vecteurs de canal h Z J en appliquant les etapes A.1, A.2 et A.3 afin d'obtenir 
les trains de symboles {b m _ k) } associe aux vecteurs de canal h sj = h r (* y ) . 
Etape n°IH.c.4 : Construction du vecteur w de I'equation (32) a partir des 

10 Etape n°lll.c.5 : Maximisation du critere Porteuse(f 0 ) de ['equation (33) pour 
obtenir f 0 . 

Etape n°m.c.6 : Application de I'equation (30) pour deduire les symboles {a m } 
des symboles {b m }. 

Etape n°lll.c.7 : Choix parmi les trains de symboles celui qui est associe au 
15 vecteur h 2 (i) de plus fort module : {a m .j}. 

Etape n°Ul.c.8 : Etapes ILb.8 jusqu'a ll.b.10 prec^demment decrites. 
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REVENDICATIONS 



10 



1 - Procede de demodulation aveugle d'une source ou emetteur de forme 
d'onde lineaire dans un systeme comportant une ou plusieurs sources et un 
reseau de capteurs et un canal de propagation caracterise en ce qu'il 
comporte au moins les etapes suivantes : 

• determiner le temps symbole T et on echantillonne a Te tel que T=ITe 
(I entier), 

• a partir des observations x(kTe), construire une observation spatio- 
temporelle z(t) dont les sources melangees sont des trains de 
symbole de l'<§metteur , 

• appliquer une methode de type ICA sur le vecteur d'observation z(t) 
pour estimer les U trains de symboies { a„* } associes aux vecteurs 

15 de canal h ..., = £..(*,). 

• ordonner les U sorties (a*,,, h ZJ ) dans le meme ordre que les entrees 
(a m -i , h z (i)) afin d'obtenir les vecteurs de canal de propagation 
h.v=A(*,) 

• determiner la phase ow associee aux sorties. 
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2 - ProcedS selon la revendication 1 caracterise en ce que I'on estime les 
parametres du canal de propagation pour determiner la frequence porteuse 
afin de compenser les trains de symboies pour les obtenir en bande de base. 

25 3 - Procede selon la revendication 1 caracterise en ce quMI comporte une 
etape d'estimation des parametres du canal de propagation en angle e p et 
retard t p .. 
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